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ABSTRACT 

We have analyzed HST/NICM0S2 FllOW, F160W, F165M, and F207M band 
images covering the central I'xl' region of the cluster associated with Mon R2 
in order to constrain the Initial Mass Function (IMF) down to 20 Mjup. The flux 
ratio between the F165M and F160W bands was used to measure the strength 
of the water band absorption feature and select a sample of 12 out of the total 
sample of 181 objects that have effective temperatures between 2700 K and 3300 
K. These objects are placed in the HR diagram together with sources observed 
by Carpenter et al. (1997) to estimate an age of ~1 Myr for the low mass cluster 
population. By constructing extinction limited samples, we are able to constrain 
the IMF and the fraction of stars with a circumstellar disk in a sample that is 90% 
complete for both high and low mass objects. For stars with estimated masses 
between 0.1 Mq and 1.0 M© for a 1 Myr population with Ay < 19 mag, we find 
that 27 ± 9% have a near-infrared excess indicative of a circumstellar disk. The 
derived fraction is similar to, or slightly lower than, the fraction found in other 
star forming regions of comparable age. We constrain the number of stars in the 
mass interval 0.08-1.0 M© to the number of objects in the mass interval 0.02- 
0.08 M© by forming the ratio, i?**=N(0.08-l M0)/N(O.O2-O.O8 Mq) for objects 
in an extinction limited sample complete for Ay < 7 mag. The ratio is found to 
be R** = 2.2 ± 1.3 assuming an age of 1 Myr, consistent with the similar ratio 
predicted by the system IMF proposed by Chabrier (2003). The ratio is similar 
to the ratios observed towards the Orion Nebula Cluster and IC 348 as well as 
the ratio derived in the 28 square degree survey of Taurus by Guieu et al. (2006). 
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Subject headings: stars: low-mass, brown dwarfs — stars: mass function — 
stars: pre-main sequence — stars: formation — open clusters and associations: 
individual(Mon R2) 

1. Introduction 

The shape of the IMF and whether or not it is universal have been central questions 
in astrophysics for more than 50 years. The shape of the stellar part of the IMF has been 
examined obscrvationally in some detail. Above 0.5 M0, the IMF derived for field stars is 
well characterized by a Salpeter power law, dN/dM cx M~"i,q;i = 2.35 (Kroupa 2002). At 
this characteristic mass, the slope of the power law changes to 0.2 = 1.35 and for masses 
below 0.1 M0 the IMF flattens further (Kroupa 2002; Chabricr 2003), although the slope 
is very uncertain (Allen et al. 2005). To search for variations in the IMF as a function 
of environment, one has to look at individual star-forming episodes, such as star clusters 
or associations. Young clusters in particular are useful for these studies because the more 
massive stars are still present in the cluster and the clusters have not experienced significant 
dynamical mass segregation. Recently, the IMF studies have been extended deep into the 
brown dwarf regime in several star forming regions, most notably the Orion Nebula Cluster 
(ONC, e.g. Hillenbrand & Carpenter 2000; Muench et al. 2002; Slesnick et al. 2004), IC 348 
(e.g. Najita et al. 2000; Luhman et al. 2003), and Taurus (e.g. Briceno et al. 2002; Luhman 
2004; Guieu et al. 2006). 

In order to search for variations in the sub-stellar IMF as a function of environment, 
a larger ensemble of young clusters is needed. Here, we present a study of the IMF in the 
Mon R2 cluster. The embedded cluster associated with Mon R2 was originally discovered 
by Beckwith et al. (1976) and is located at a distance of 830±50 pc (Herbst & Racine 1976). 
The most massive star in the cluster is ~ 10 M© (Carpenter et al. 1997; Massi et al. 1985). 
Carpenter et al. (1997) estimated Mon R2 to contain at least 475 stars within a 3.6x3.6 
pc region, significantly less than the ~ 3500 stars found within a 2.5 pc radius of the ONC 
(Hillenbrand & Hartmann 1998). IC 348 contains ~ 400 members including brown dwarfs 
(Luhman et al. 2003). Thus, Mon R2 has a richness intermediate between the two well 
studied clusters. The molecular cloud has been shown by Choi et al. (2000) to have a central 
density of rif. ~ lO^cm^"^, much higher than typical densities in the Taurus molecular cloud 
cores {ric ~ lO^cm^^, see e.g. Onishi et al. (2002)). If the density of the molecular cores 
have an impact on the IMF as suggested by e.g. Bate & Bonnell (2005) and Goodwin et 
al. (2004), we would expect the IMF in Mon R2 to have relatively more brown dwarfs than 
Taurus. 
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Previously, Carpenter et al. (1997) imaged Mon R2 down to the brown dwarf limit and 
obtained spectra for a subset of the stars. They estabhshed the age of the chister to be < 
3 Myr, and that the majority of the low mass stars clustered around the 1 Myr isochrone. 
Adopting an age of 1 Myr for the stellar population, they found the ratio of high mass to 
low mass stars N(l-10 M0)/N(O.1-1 M0)= 0.1, which is consistent with other nearby star 
forming regions (Meyer et al. 2000) and with a Miller & Scalo (1979) field star IMF. Despite 
the fact that Mon R2 is relatively nearby, the brown dwarf content has not been investigated 
in a systematic manner. Here, we present HST/NICMOS 2 near-infrared imaging of the 
central 0.24pcx0.24pc (1' x 1') of the cluster. Our main goal is to constrain the low-mass 
stellar and brown dwarf content in the cluster and compare the results to those derived for 
other young clusters and the field. In addition, we estimate the fraction of stars with disks 
within our surveyed region. 

The paper is structured as follows: In Section 2, we present the observations, the data 
reduction and the photometry. Section 3 presents the results from the color-magnitude 
and color-color diagrams for the cluster. We place objects in the HR diagram based on the 
estimate of the effective temperature derived from the F160W and F165M band observations. 
In Section 4, we discuss our constraints on the IMF and disk frequency based on extinction 
limited samples constructed from our observations. Section 5 compares our derived ratios 
with the similar ratios for other star forming regions. Finally, we conclude in Section 6 with 
a summary of our results. 

2. Observations and reduction 
2.1. The dataset 

HST/NICM0S2 observations of the central I'xl' of Mon R2 were obtained in Cycle 
7 under program number 7417. The field of view for NICMOS 2 is 19'.'2 square and the 
pixel scale is 0'.'075. A 4 x 4 grid was observed in the FllOW, F160W, F165M, and F207M 
filters, where each position in the grid was observed twice in each filter with a dither offset 
of ~ 20 pixels in both RA and DEC. The total integration time per position in the grid is 
512 seconds for FllOW, 352 seconds for F160W, 576 seconds for F165M, and 288 seconds 
for F207M. All the images were obtained in non-destructive readout mode with a total of 
16 readouts in each exposure. The full width at half maxima of the point spread function 
are 0'/12, 0'/15, and 0'.'20, for the FllOW, F160W, and F207M bands respectively. Table 1 
contains the observing log. 

The images were processed using a combination of IRAF scripts, C programs, and IDL 
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scripts. The basic reduction was performed using the NICRED procedure (Lehar et al. 
2000) which performs hnearity corrections, cosmic ray rejection, determines photon rates, 
does dark subtraction, and fiat fielding. Synthetic darks were used based on the software 
NICRED. The pedestal effect was removed from each frame also using NICRED. Bad pixels 
were identified and corrected using the IRAF task f ixpix. Finally, for each filter, the images 
were registered and a mosaic was created using the average of the two dithered exposures 
at each position in the mosaic. A color composite of the observed region is shown in Fig. 1, 
where FllOW is blue, F160W is green, and F207M red. 



Table 1: Observing log for the NICM0S2 observations. 

Filter Obs. date Exp. time 
FllOW Dec 1 1997 2x256 
F160W Dec 1 1997 2x176 
F165M Dec 5 1997 2x288 
F207M Dec 5 1997 2x144 
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Fig. 1.— HST/NICMOS 2 FllOW (blue), F160W (green), and F207M (red) color- 
composite of the central Ixl'of Mon R2. The faintest (red) stars seen have a magnitude of 
fnF2Q7M ~19.5 mag. North is up and east left. The coordinates of the bright red source in 
the center of the image is (RA,DEC)=(06:07:45.77,-06:22:53.6) 
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2.2. Source detection and photometry 

Point sources in each of the mosaics were identified using the IRAF implementation 
of daof ind with a 5a detection threshold. False detections in the form of noise peaks or 
structures in the point spread function were eliminated by comparing the source lists from 
each filter and a visual inspection of each source detected by daof ind. Since the location of 
features in the PSF for a given star is dependent on wavelength, the use of different filters 
helps in separating faint stars next to bright stars from structure in the bright star PSF. 

Aperture photometry was performed using the apphot package in IRAF. An aperture 
of 3.5 pixel radius was used to measure the fiux of the object and the background was 
measured in an annulus from 10 to 15 pixels. For 10 objects an annulus of 2 pixel was used 
to measure the object flux because they were close to bright stars. The choice of a relatively 
large aperture for the majority of the sources was to limit photometric errors due to the 
variations in the PSF across the field of view of the NIC2 camera. Some 40 stars with no 
(or very little) nebulosity associated with them were used to determine aperture corrections 
for the photometry in each filter. The NICMOS photometry was measured relative to a 
10 pixel radius aperture and is presented in the VEGAMAG system. The astrometry was 
done relative to Carpenter et al. (1997). Some 25 stars in common were used to fit a 
second order polynomial in both RA and DEC. The final RMS was O'.'l. A comparison with 
2MASS showed an offset by 0'.'5 in RA and -0'.'3 in DEC relative to our astrometry. We 
corrected the astrometry for this offset and the RMS compared to 2MASS is then 0'.'2 for 
8 relatively bright isolated sources. Table 2 presents the photometry and coordinates for 
the 181 sources detected together with 2MASS IDs where available. Where more than one 
source was within an arcsecond in the HST/NICMOS survey, we have chosen the brightest 
source in our observations as the 2MASS source. In four cases, two sources of roughly equal 
brightness in the NICMOS data were identified as one object by 2MASS. Both objects have 
been given the 2MASS ID and are marked with an *. 



Table 2. The HST/NICMOS photometry for Mon R2. 



ID RA (2000) DEC (2000) FllOW F160W F165M F207M 2MASS ID 



1 


6:07:44.98 


-6:23; 


;25.8 


23, 


,11± 


0.51 


17.32 


± 


0.05 


17.10 


±0.05 


14.47 


± 


0.06 


06074492-0623257 


2 


6:07:45.63 


-6:23: 


:25.1 


24, 


,55± 


1.21 


19.77 


± 


0.07 


19.41 


±0.07 


17.07 


± 


0.07 




3 


6:07:43.89 


-6:23: 


:26.8 


18, 


,58± 


0.03 


16.20 


± 


0.04 


15.90 


±0.04 


15.48 


± 


0.06 




4 


6:07:44.82 


-6:23: 


:25.8 


17, 


,34± 


0.03 


14.77 


± 


0.03 


14.55 


±0.03 


13.89 


± 


0.05 


06074483-0623260 


5 


6:07:44.26 


-6:23; 


;25.5 


16, 


.12± 


0.04 


13.81 


± 


0.04 


13.58 


±0.04 


13.04 


± 


0.06 


06074426-0623256 


6 


6:07:47.67 


-6:23; 


;22.0 


20, 


,17± 


0.04 


18.45 


± 


0.04 


18.15 


±0.04 


17.90 


± 


0.08 




7 


6:07:46.08 


-6:23; 


;21.8 


18, 


.62± 


0.03 


15.92 


± 


0.03 


15.68 


±0.03 


15.01 


± 


0.05 




8 


6:07:44.98 


-6:23; 


:21.9 


21, 


.13± 


0.05 


17.56 


± 


0.03 


17.29 


±0.03 


16.17 


± 


0.05 




9 


6:07:47.23 


-6:23; 


18.8 


20, 


,88± 


0.05 


18.74 


± 


0.04 


18.51 


±0.04 


18.03 


± 


0.08 




10 


6:07:44.53 


-6:23; 


;21.4 


17, 


.75± 


0.04 


14.96 


± 


0.04 


14.67 


±0.04 


13.63 


± 


0.06 


06074454-06232 14 


11 


6:07:45.31 


-6:23; 


;20.2 


15, 


,51± 


0.04 


12.99 


± 


0.04 


12.76 


±0.04 


12.13 


± 


0.06 


06074531-0623202 


12 


6:07:43.82 


-6:23; 


;21.5 


21, 


,76± 


0.11 


18.46 


± 


0.04 


18.14 


±0.04 


17.08 


± 


0.07 




13 


6:07:45.21 


-6:23; 


;19.5 


18, 


.62± 


0.04 


16.35 


± 


0.04 


16.06 


±0.04 


15.62 


± 


0.07 




14 


6:07:46.64 


-6:23; 


:17.3 


15, 


,26± 


0.04 


12.80 


± 


0.04 


12.90 


±0.04 


12.09 


± 


0.06 


06074665-0623175 


15 


6:07:45.68 


-6:23; 


:17.7 


18, 


,91± 


0.04 


16.26 


± 


0.04 


16.02 


±0.04 


15.37 


± 


0.07 




16 


6:07:43.96 


-6:23; 


18.5 


17, 


,00± 


0.04 


14.52 


± 


0.04 


14.29 


±0.04 


13.71 


± 


0.06 


06074396-0623184 


17 


6:07:44.10 


-6:23; 


;18.2 


21, 


,25± 


0.06 


18.47 


± 


0.04 


18.28 


±0.04 


17.26 


± 


0.07 




18 


6:07:45.49 


-6:23; 


;16.0 


16, 


,32± 


0.04 


14.06 


± 


0.04 


13.84 


±0.04 


13.43 


± 


0.06 


06074549-0623160 


19 


6:07:44.95 


-6:23; 


;16.4 


22, 


,56± 


0.16 


19.34 


± 


0.05 


19.06 


±0.05 


17.05 


± 


0.05 




20 


6:07:46.31 


-6:23; 


;13.7 


18, 


.43± 


0.04 


16.00 


± 


0.04 


15.72 


±0.04 


15.21 


± 


0.07 


06074631-0623128 


21 


6:07:44.62 


-6:23; 


:15.3 


22, 


,75± 


0.15 


17.75 


± 


0.04 


17.42 


±0.04 


15.56 


± 


0.06 




22 


6:07:45.65 


-6:23; 


:14.1 


18, 


.46± 


0.04 


16.29 


± 


0.04 


16.06 


±0.04 


15.71 


± 


0.07 




23 


6:07:44.37 


-6:23; 


15.4 








20.90 


± 


0.17 


20.79 


±0.18 


18.54 


± 


0.12 




24 


6:07:44.48 


-6:23; 


;15.2 








20.43 


± 


0.12 


20.40 


±0.14 


17.49 


± 


0.08 




25 


6:07:47.74 


-6:23: 


;11.8 


16, 


,98± 


0.04 


15.01 


± 


0.04 


14.77 


±0.04 


14.45 


± 


0.06 




26 


6:07:47.12 


-6:23; 


;12.2 


21, 


.33± 


0.07 


18.75 


± 


0.05 


18.47 


±0.05 


17.74 


± 


0.09 




27 


6:07:46.69 


-6:23; 


:12.6 


15, 


.84± 


0.04 


13.58 


± 


0.04 


13.35 


±0.04 


12.87 


± 


0.06 


06074669-0623125 


28 


6:07:44.80 


-6:23; 


:14.4 


16, 


,95± 


0.03 


13.79 


± 


0.03 


13.59 


±0.03 


12.59 


± 


0.05 


06074482-0623144 


29 


6:07:46.33 


-6:23; 


;12.1 


20, 


.55± 


0.05 


16.52 


± 


0.04 


16.21 


±0.04 


14.37 


± 


0.06 




30 


6:07:45.46 


-6:23; 


;12.9 


21, 


.62± 


0.14 


18.59 


± 


0.08 


18.33 


±0.08 


17.50 


± 


0.12 




31 


6:07:47.26 


-6:23; 


;10.5 


24, 


,16± 


0.67 


21.84 


± 


0.69 


22.11 


±1.15 


18.99 


± 


0.24 




32 


6:07:43.95 


-6:23; 


;13.0 


21, 


.38± 


0.07 


16.62 


± 


0.04 


16.40 


±0.04 


14.10 


± 


0.06 


06074393-0623130 


33 


6:07:46.23 


-6:23; 


:10.4 


















17.91 


± 


0.18 




34 


6:07:45.73 


-6:23; 


:10.9 


20, 


,83± 


0.05 


17.38 


± 


0.03 


17.06 


±0.03 


15.79 


± 


0.05 


06074579-0623113 


35 


6:07:45.11 


-6:23; 


:10.8 


20, 


,79± 


0.06 


16.52 


± 


0.04 


16.07 


±0.04 


14.16 


± 


0.06 


06074593-0623110 


36 


6:07:45.80 


-6:23; 


;09.6 


20, 


.79± 


0.04 


16.93 


± 


0.02 


16.63 


±0.02 


15.24 


± 


0.04 


06074500-0623113 


37 


6:07:45.96 


-6:23; 


;09.2 


21, 


,37± 


0.14 


17.80 


± 


0.04 


17.63 


±0.04 


15.67 


± 


0.04 




38 


6:07:46.16 


-6:23; 


;09.0 








19.18 


± 


0.09 


19.02 


±0.09 


16.35 


± 


0.06 




39 


6:07:46.86 


-6:23; 


;07.9 


16, 


.73± 


0.02 


14.64 


± 


0.02 


14.39 


±0.02 


13.90 


± 


0.03 


06074686-0623076 


40 


6:07:43.99 


-6:23; 


:10.6 


19, 


,51± 


0.03 


16.63 


± 


0.03 


16.36 


±0.03 


15.60 


± 


0.05 




41 


6:07:43.68 


-6:23; 


10.4 


17 


,50dz 


0.03 


15.12 


± 


0.03 


14.84 


±0.03 


14.29 


± 


0.05 


06074370-0623107 


42 


6:07:45.16 


-6:23; 


:08.7 


20, 


,17± 


0.04 


16.42 




0.03 


16.16 


±0.03 


14.96 


± 


0.05 


06074516-0623103 


43 


6:07:45.42 


-6:23; 


;08.2 


19, 


,23± 


0.03 


15.62 


± 


0.03 


15.11 


±0.03 


13.86 


± 


0.05 


06074542-0623092 


44 


6:07:44.81 


-6:23; 


;08.4 


24, 


,04± 


0.68 


19.33 


± 


0.04 


19.08 


±0.03 


16.69 


± 


0.04 




45 


6:07:45.89 


-6:23; 


;07.3 


22, 


,60± 


0.16 


19.47 


± 


0.04 


19.21 


±0.04 


17.65 


± 


0.05 




46 


6:07:46.64 


-6:23; 


;05.9 


20, 


,61± 


0.12 


17.02 


± 


0.07 


16.78 


±0.07 


15.41 


± 


0.08 


06074661-0623058 
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Table 2 — Continued 



ID RA (2000) DEC (2000) FllOW F160W F165M F207M 2MASS ID 



06074784-0623038 
06074503-0623060 

06074615-0623046 
06074470-0623044 

06074807-0623017 
06074470-0623044 



06074622-0623022 



06074495-0623010 



47 


6:07:47.25 


-6; 


:23: 


:05. 


,1 


21, 


.95± 


0, 


,15 


19, 


,59 


± 





.08 


19, 


,27 


±0, 


,08 


18, 


,18 


± 


0.15 


48 


6:07:46.62 


-6: 


:23: 


:05, 


,7 


21, 


,01± 


0, 


,19 


18, 


,23 


± 


0, 


,09 


17 


,94 


±0, 


,09 


16, 


,45 


± 


0.08 


49 


6:07:45.64 


-6: 


:23: 


:06, 


,5 


20, 


,92± 


0, 


,05 


16, 


,55 


± 


0, 


,03 


16, 


,23 


±0, 


,03 


14, 


,78 


± 


0.05 


50 


6:07:47.84 


-6; 


:23: 


:03, 


,7 


16, 


,52± 


0, 


,03 


13, 


,82 




0, 


,03 


13, 


,58 


±0, 


,03 


12, 


,82 


± 


0.05 


51 


6:07:45.03 


-6: 


;23: 


;06. 


,5 


20, 


.49± 


0, 


,04 


16, 


,78 


± 


0, 


.03 


16, 


,68 


±0, 


,03 


15, 


,49 


± 


0.05 


52 


6:07:45.08 


-6; 


:23: 


:06. 


.0 


21, 


,60± 


0, 


,11 


17, 


,03 


± 


0, 


.04 


16, 


,71 


±0, 


,04 


15, 


,18 


± 


0.07 


53 


6:07:46.13 


-6; 


:23: 


:04. 


.5 


17, 


,19± 


0, 


,04 


12, 


,83 


± 


0, 


.04 


12, 


,57 


±0, 


,04 


10, 


,78 


± 


0.06 


54 


6:07:44.62 


-6; 


:23: 


:05. 


,3 


19, 


.96± 


0, 


,04 


16, 


,66 


± 





.04 


16, 


,39 


±0, 


,04 


15, 


,25 


± 


0.07 


55 


6:07:48.07 


-6: 


:23: 


:01, 


,9 


23, 


,92± 


0, 


,81 


18, 


,65 


± 


0, 


,04 


18, 


,34 


±0, 


,04 


16, 


,13 


± 


0.05 


56 


6:07:44.80 


-6: 


:23: 


:05, 


1 


17 


,33± 


0, 


,04 


15, 


,51 


± 


0, 


.04 


15, 


,27 


±0, 


,04 


15, 


,19 


± 


0.07 


57 


6:07:45.88 


-6: 


:23: 


:03, 


,5 


22, 


,96± 


0, 


,38 


18, 


,92 


± 


0, 


.04 


18, 


,56 


±0, 


,04 


16, 


,06 


± 


0.05 


58 


6:07:47.67 


-6: 


;23: 


;01. 


,1 


18, 


.91± 


0, 


,04 


16, 


,95 


± 





.07 


16, 


,64 


±0, 


,07 


16, 


,56 


± 


0.28 


59 


6:07:45.26 


-6: 


:23: 


:03. 


.2 


22, 


,02± 


0, 


,10 


17, 


,65 


± 


0, 


.04 


17, 


,58 


±0, 


,04 


15, 


,31 


± 


0.06 


60 


6:07:44.62 


-6: 


:23: 


:03. 


.9 










19, 


,66 


± 


0, 


.14 


19, 


,21 


±0, 


,11 


16, 


,52 


± 


0.08 


61 


6:07:46.22 


-6; 


:23: 


:01. 


,9 


12, 


.34± 


0, 


,04 


11, 


,57 


± 





.04 


11, 


,45 


±0, 


,04 


11, 


,63 


± 


0.06 


62 


6:07:46.02 


-6: 


:23: 


:01, 


,6 










18, 


,52 


± 


0, 


.07 


18, 


,31 


±0, 


,06 


14, 


,78 


± 


0.06 


63 


6:07:45.67 


-6: 


:23: 


:01, 


,9 


23, 


,57± 


0, 


,47 


19, 


,77 


± 


0, 


.09 


19, 


,48 


±0, 


,09 


17 


,01 


± 


0.08 


64 


6:07:45.85 


-6: 


:23: 


:01, 


,4 


25, 


,11± 


3, 


,73 


21, 


,61 


± 


0, 


.29 


20, 


,85 


±0, 


,18 


18, 


,57 


± 


0.13 


65 


6:07:45.99 


-6: 


;23: 


;01. 


,1 


21, 


.89± 


0, 


,17 


20, 


,78 


± 





.19 


20, 


,26 


±0, 


,17 


18, 


,08 


± 


0.11 


66 


6:07:45.04 


-6: 


:23: 


:01. 


,8 


21, 


,59± 


0, 


,18 


18, 


,76 


± 





.11 


18, 


,43 


±0, 


,09 


17, 


,49 


± 


0.13 


67 


6:07:46.16 


-6; 


;23: 


;00. 


,4 


22, 


.44± 


0, 


,38 


18, 


,79 


± 





.07 


18, 


,41 


±0, 


,06 


16, 


,12 


± 


0.07 


68 


6:07:45.65 


-6: 


:23: 


:00, 


,8 


26, 


.05± 


3, 


,84 


21, 


,88 


± 





.56 


21, 


,53 


±0, 


,46 


18, 


,51 


± 


0.17 


69 


6:07:46.01 


-6: 


:23: 


:00, 


,3 


23, 


,67± 


1, 


,21 












23, 


,56 


±4, 


,82 


21, 


,81 


± 


3.48 


70 


6:07:45.98 


-6: 


;22: 


;59. 


,9 


23, 


.32± 


0, 


,82 




















19, 


,46 


± 


0.46 


71 


6:07:45.71 


-6: 


:23: 


:00. 


,0 


22, 


,59± 


0, 


,22 


18, 


,03 


± 


0, 


.03 


17, 


,69 


±0, 


,03 


15, 


,05 


± 


0.05 


72 


6:07:45.64 


-6: 


:22: 


:59. 


,4 


23, 


,07± 


0, 


,23 


21, 


,15 


± 





.29 


20, 


,38 


±0, 


,16 


19, 


,48 


± 


0.43 


73 


6:07:46.75 


-6; 


:22: 


:57. 


,6 


21, 


.64± 


0, 


,17 


18, 


,45 


± 





.10 


18, 


,07 


±0, 


,10 


16, 


,49 


± 


0.10 


74 


6:07:46.72 


-6: 


:22: 


:57, 


,5 


20, 


,21± 


0, 


,08 


16, 


,54 


± 


0, 


,07 


16, 


,19 


±0, 


,07 


14, 


,95 


± 


0.07 


75 


6:07:47.53 


-6: 


:22: 


:56, 


,5 


17 


,33± 


0, 


,04 


15, 


,11 


± 


0, 


.04 


14, 


,84 


±0, 


,04 


14, 


,35 


± 


0.07 


76 


6:07:46.89 


-6: 


:22: 


;57. 


,0 


20, 


.58± 


0, 


,06 


18, 


,54 


± 


0, 


.04 


18, 


,25 


±0, 


,04 


17 


,72 


± 


0.08 


77 


6:07:47.84 


-6: 


;22: 


:56. 


,1 


14, 


.23± 


0, 


,04 


10, 


,43 


± 





.04 


10, 


,16 


±0, 


,04 


8, 


,21 


± 


0.06 


78 


6:07:47.74 


-6; 


:22: 


:56. 


.1 


17, 


,89± 


0, 


,05 


14, 


,75 


± 


0, 


.04 


14, 


,45 


±0, 


,04 


13, 


,04 


± 


0.06 


79 


6:07:46.28 


-6: 


:22: 


:57. 


.1 










18, 


,05 


± 





.04 


18, 


,22 


±0, 


,05 


15, 


,31 


± 


0.06 


80 


6:07:47.89 


-6; 


:22: 


:55. 


,2 


15, 


.25± 


0, 


,04 


11, 


,91 


± 





.04 


11, 


,68 


±0, 


,04 


9, 


,73 


± 


0.06 


81 


6:07:46.43 


-6: 


:22: 


:56, 


,3 


22, 


,72± 


0, 


,25 


19, 


,58 


± 





,07 


19, 


,19 


±0, 


,06 


17 


,69 


± 


0.08 


82 


6:07:47.68 


-6: 


:22: 


:55, 


,0 


15, 


,70± 


0, 


,04 


13, 


,66 


± 


0, 


,04 


13, 


,34 


±0, 


,04 


12, 


,78 


± 


0.06 


83 


6:07:46.89 


-6: 


:22: 


:55, 


,6 


19, 


,27± 


0, 


,04 


18, 


,27 


± 


0, 


.04 


18, 


,07 


±0, 


,04 


18, 


,15 


± 


0.08 


84 


6:07:45.02 


-6: 


;22: 


;56. 


,9 


22, 


.39± 


0, 


,12 


18, 


,05 


± 





.05 


17 


,76 


±0, 


,05 


16, 


,14 


± 


0.07 


85 


6:07:47.89 


-6: 


:22: 


:54. 


,0 


17, 


,08± 


0, 


,04 


13, 


,85 


± 





.04 


13, 


,71 


±0, 


,04 


12, 


,42 


± 


0.08 


86 


6:07:46.69 


-6: 


:22: 


:54. 


.8 


19, 


,66± 


0, 


,05 


15, 


,88 


± 


0, 


.04 


15, 


,61 


±0, 


,04 


14, 


,42 


± 


0.06 


87 


6:07:47.38 


-6; 


:22: 


:53. 


,9 


17, 


.61± 


0, 


,04 


15, 


,05 


± 





.04 


14, 


,81 


±0, 


,04 


14, 


,03 


± 


0.06 


88 


6:07:45.45 


-6: 


;22: 


;55, 


,7 


23, 


,27± 


0, 


,34 


19, 


,46 







,09 


19, 


,08 


±0, 


,08 


17 


,65 


± 


0.11 


89 


6:07:45.72 


-6; 


;22: 


:55. 


,0 


23, 


,12± 


0, 


,30 


18, 


,90 


± 


0, 


.04 


18, 


,62 


±0, 


,04 


15, 


,44 


± 


0.05 


90 


6:07:47.67 


-6; 


:22: 


:53. 


.0 


20, 


.85± 


0, 


,08 


18, 


,02 


± 


0, 


.10 


17, 


,74 


±0, 


,10 


16, 


,04 


± 


0.09 


91 


6:07:46.03 


-6: 


:22: 


:54. 


,5 


25, 


.26± 


2, 


,81 


19, 


,39 


± 





.09 


19, 


,00 


±0, 


,08 


16, 


,60 


± 


0.07 


92 


6:07:46.37 


-6: 


:22: 


:54. 


,0 


22, 


.33± 


0, 


,20 


19, 


,59 


± 





.09 


19, 


,24 


±0, 


,08 


17, 


,83 


± 


0.09 



06074674-0622562 



06074786-0622559 



06074488-0622563 



-9- 



Table 2 — Continued 



ID RA (2000) DEC (2000) 



FllOW 



F160W 



F165M 



F207M 



2MASS ID 



93 


6 


07:46 


10 


-6 


22 


54 


1 










18 


57 


± 





04 


18 


40 


±0 


04 


15 


39 


± 





05 


94 


6 


07:45 


64 


-6 


22 


54 


3 


24 


30± 





67 


18 


46 


± 





04 


18 


30 


±0 


04 


15 


97 


± 





05 


95 


6 


07:44 


14 


-6 


22 


55 


5 


23 


41± 





29 


20 


73 


± 





17 


20 


32 


±0 


15 


19 


14 


± 





24 


96 


6 


07:45 


80 


-6 


22 


53 


3 


19 


54± 





02 


14 


59 


± 





02 


14 


52 


±0 


02 


9 


77 


± 





03 


97 


6 


07:46 


01 


-6 


22 


53 





23 


28± 





36 


19 


93 


± 





09 


19 


75 


±0 


08 


16 


99 


± 





06 


98 


6 


07:46 


60 


-6 


22 


52 


4 


19 


67± 





03 


16 


12 


± 





03 


15 


77 


±0 


03 


14 


69 


± 





05 


99 


6 


07:45 


69 


-6 


22 


52 


4 


22 


92± 





31 


18 


80 


± 





04 


18 


87 


±0 


04 


15 


42 


± 





04 


100 


6 


07:45 


75 


-6 


22 


51 


7 


22 


53± 





19 


20 


97 


± 





19 


21 


18 


±0 


34 


16 


68 


± 





07 


101 


6 


07:46 


42 


-6 


22 


50 


8 


22 


43± 





14 


19 


91 


± 





13 


19 


75 


d=0 


13 


18 


33 


± 





19 


102 


6 


07:44 


54 


-6 


22 


52 


6 


23 


87± 





60 


20 


54 


± 





13 


20 


20 


±0 


12 


17 


76 


± 





06 


103 


6 


07:46 


47 


-6 


22 


50 


4 


21 


84± 





08 


17 


90 


± 





03 


17 


60 


±0 


03 


16 


03 


± 





05 


104 


6 


07:43 


64 


-6 


22 


52 


5 


16 


40± 





04 


13 


84 


± 





04 


13 


64 


±0 


04 


12 


80 


± 





06 


105 


6 


07:46 


38 


-6 


22 


49 


5 


21 


53± 





08 


18 


22 


± 





04 


17 


83 


±0 


04 


16 


50 


± 





06 


106 


6 


07:46 


18 


-6 


22 


49 


6 


23 


01± 





24 


19 


87 


± 





16 


19 


50 


±0 


14 


17 


52 


± 





16 


107 


6 


07:45 


60 


-6 


22 


50 


1 


21 


99± 





10 


17 


74 


± 





03 


17 


48 


±0 


03 


15 


70 


± 





05 


108 


6 


07:46 


26 


-6 


22 


49 


4 


23 


11± 





28 


17 


94 


± 





05 


17 


61 


±0 


05 


15 


46 


± 





07 


109 


6 


07:45 


70 


-6 


22 


49 


4 


21 


68± 





09 


17 


13 


± 





03 


16 


79 


±0 


03 


15 


11 


± 





05 


110 


6 


07:44 


10 


-6 


22 


50 


9 


22 


68± 





25 


19 


63 


± 





07 


19 


25 


±0 


06 


18 


04 


± 





09 


111 


6 


07:43 


60 


-6 


22 


51 


3 


15 


61± 





04 


13 


15 


± 





04 


12 


93 


±0 


04 


12 


43 


± 





06 


112 


6 


07:45 


56 


-6 


22 


49 


3 


22 


61± 





19 


17 


51 


± 





04 


17 


19 


±0 


04 


14 


94 


± 





07 


113 


6 


07:46 


33 


-6 


22 


46 


8 


20 


40± 





05 


17 


12 


± 





04 


16 


78 


±0 


04 


15 


87 


± 





07 


114 


6 


07:45 


30 


-6 


22 


47 


1 


22 


72± 





20 


18 


70 


± 





06 


18 


41 


±0 


05 


16 


19 


± 





07 


115 


6 


07:46 


27 


-6 


22 


45 


6 


21 


50± 





11 


17 


99 


± 





06 


17 


72 


±0 


05 


16 


40 


± 





08 


116 


6 


07:46 


49 


-6 


22 


45 





22 


63± 





17 


17 


62 


± 





04 


17 


41 


±0 


04 


15 


09 


± 





07 


117 


6 


07:45 


07 


-6 


22 


46 


3 


20 


86± 





06 


16 


71 


± 





04 


16 


48 


±0 


04 


14 


87 


± 





06 


118 


6 


07:44 


97 


-6 


22 


46 


4 


21 


44± 





08 


17 


52 


± 





04 


17 


23 


±0 


04 


15 


92 


± 





07 


119 


6 


07:45 


16 


-6 


22 


45 


7 


22 


76± 





23 


20 


80 


± 





16 


20 


37 


±0 


13 


19 


45 


± 





34 


120 


6 


07:43 


74 


-6 


22 


46 


7 


16 


92± 





09 


14 


71 


± 





09 


14 


40 


±0 


09 


14 


00 


d= 





10 


121 


6 


07:43 


75 


-6 


22 


46 


4 


17 


77± 





09 


15 


41 


± 





09 


15 


15 


±0 


09 


14 


70 


± 





10 


122 


6 


07:44 


45 


-6 


22 


44 


8 


22 


07± 





13 


17 


42 


± 





04 


17 


14 


±0 


04 


15 


49 


± 





06 


123 


6 


07:46 


40 


-6 


22 


42 


7 


17 


74± 





04 


15 


21 


± 





04 


15 


00 


±0 


04 


14 


42 


± 





06 


124 


6 


07:43 


55 


-6 


22 


44 


6 


19 


28± 





09 


16 


04 


± 





09 


15 


74 


±0 


09 


14 


58 


± 





10 


125 


6 


07:43 


56 


-6 


22 


44 


6 


18 


74± 





09 


15 


46 


± 





09 


15 


17 


±0 


09 


14 


20 


± 





10 


126 


6 


07:43 


90 


-6 


22 


44 


2 


22 


41± 





20 


20 


15 


± 





07 


19 


84 


±0 


07 


18 


93 


± 





12 


127 


6 


07:46 


94 


-6 


22 


40 


9 


21 


87± 





09 


17 


79 


± 





03 


17 


49 


±0 


03 


16 


10 


± 





05 


128 


6 


07:46 


27 


-6 


22 


41 





23 


00± 





27 


20 


57 


± 





20 


20 


00 


±0 


16 


18 


91 


± 





33 


129 


6 


07:44 


99 


-6 


22 


42 


2 


24 


44± 


1 


09 


19 


35 


± 





07 


18 


97 


±0 


06 


17 


48 


± 





08 


130 


6 


07:45 


81 


-6 


22 


41 


4 


24 


99± 


1 


97 


19 


06 


± 





05 


18 


88 


±0 


05 


16 


53 


± 





05 


131 


6 


07:44 


76 


-6 


22 


41 


5 


24 


28± 





96 


20 


88 


± 





12 


20 


82 


±0 


13 


18 


33 


± 





08 


132 


6 


07:47 


43 


-6 


22 


38 


5 


24 


01± 





44 


23 


40 


± 


1 


12 


24 


70 


±4 


58 


18 


28 


± 





09 


133 


6 


07:44 


48 


-6 


22 


41 


2 


18 


30± 





04 


15 


11 


± 





04 


14 


88 


±0 


04 


13 


93 


± 





06 


134 


6 


07:46 


36 


-6 


22 


38 


8 


21 


17± 





07 


17 


37 


± 





04 


17 


03 


±0 


04 


15 


79 


± 





07 


135 


6 


07:44 


96 


-6 


22 


40 





16 


57± 





03 


13 


86 


± 





03 


13 


90 


±0 


03 


13 


02 


± 





05 


136 


6 


07:46 


32 


-6 


22 


38 


4 


18 


28± 





02 


15 


13 


± 





02 


14 


82 


±0 


02 


13 


43 


± 





03 


137 


6 


07:46 


59 


-6 


22 


37 


2 


22 


06± 





12 


19 


68 


± 





10 


19 


47 


±0 


11 


18 


05 


± 





12 


138 


6 


07:43 


70 


-6 


22 


40 


1 


20 


22± 





04 


15 


90 


± 





03 


15 


65 


±0 


03 


13 


46 


± 





05 



06074443-0622523 
06074361-0622519* 



06074361-0622519* 
06074636-0622467 



06074374-0622467 



06074640-0622432 
06074356-0622448* 
06074356-0622448* 



06074449-0622414 
06074496-0622401 

06074371-0622402 



Table 2 — Continued 



ID RA (2000) DEC (2000) FllOW F160W F165M F207M 2MASS ID 



139 


6 


07:45.78 


-6:22 


37.6 


17.95ib 


Q 


03 


15.52 


-|- 


0.03 


15.06 


±0.03 


14.51 


-|- 


0.05 




140 


6 




-6:22 


37.7 


22.41it 


Q 


1^ 


19.36 


-(- 


0.05 


19.13 


±0.05 


17.40 


-|- 


0.07 




141 


fi 
\j 


n7'Afi 


-6*22 


00.0 


22 16ih 


n 
u 


14 




-1- 


u.uo 


19.76 


±0.08 




-1- 


u.uo 


UOU / f±000-UO^ZO / u 


142 


Q 


U 1 .^u.uu 


-6:22 


37.4 


17.58± 





03 


13.26 


-|- 


0.03 


13.04 


±0.03 


11.03 


-|- 


0.05 


nfin7A.'^fi'?-nfi99^><n 

UUU 1 1UUO uu^^oou 


143 


Q 


07:47.01 


-6:22 


35.4 


19.83it 





03 


16.74 


-|- 


0.03 


16.42 


±0.03 


15.62 


-|- 


0.05 




144 


6 


u 1 .10 . 0'-t 


-6:22 


38.4 


22.91i 


Q 


22 


19.27 


-j- 


0.03 


18.79 


±0.03 


17.01 


-|- 


0.05 




145 


6 


07:47.48 


-6:22 


34.4 


22.98it 





17 


20.56 




0.08 


20.13 


±0.06 


19.66 


-|- 


0.25 




146 


6 


07:45.74 


-6:22 


35.3 


22.73i 


Q 


18 


18.24 




0.02 


18.03 


±0.02 


15.56 


-|- 


0.03 


nfinvzi^fi^ nfi99'?^^ 


147 


6 


07-4'^ ^'^ 
u 1 .to.ou 


-6:22 


37.0 


18.23it 


Q 


04 


15.47 


-(- 


0.04 


15.11 


±0.04 


14.51 


-|- 


0.06 


nfin7zL'?^^-nfi99'?7i 

UUU 1 'iOOU uu^^o 1 i 


148 


o 




-6*22 


00.0 


1 « OQ-I- 


n 
u 




±0.00 


-1- 




13.41 


±0.04 


12 45 


-L 

in 


u.uu 


nfin7Afifin nfi99'?'?f^ 

UUU / lOUU-UO^ZOOO 


1 


o 




-6"22 


00. u 










-L 

in 




15.92 


±0.04 




-1- 


u.uu 




150 


6 


07:47.33 


-6:22 


32.0 








22.20 


-j- 


0.44 


21.99 
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The photometry has been converted into the CIT system using 41 stars in common with 
Carpenter et al. (1997). The stars covered a color range 0.4 < mpiQQiv — mp207M < 3.2, 
where m^igow^ and mp2mM are the magnitudes in the F160W and F207M bands, respectively. 
Least-square fits were performed to convert the NICMOS photometry into the CIT system, 
resulting in fits given by: 

mj = mpiiQw - (0.789 ± 0.025) 

- (0.200 ± 0.027) X {niFimw - 207m) (1) 
rriH = mpimw - (0.089 ± 0.015) 

- (0.264 ± 0.013) X {mpimw - mpimM) (2) 
rriK = mi.207M - (0.214 ±0.018) 

- (0.306 ± 0.012) X (mpieoM - mF207M). (3) 

The photometry reaches roughly 2.5 mag deeper in the J band than the observations by 
Carpenter et al. (1997). 



2.3. Completeness limits 

Artificial star experiments have been performed in order to assess the completeness 
limits of the data. Synthetic PSFs have been created for each filter, using the Tinytim 
software (Krist et al. 1998), and a total of 160 artificial stars have been placed in each 
mosaic. The fraction of artificial stars recovered as a function of magnitude in each filter 
is presented in Table 3. All the recovery fractions arc for a 5(7 threshold and were subject 
to visual inspection in a way identical to that used to detect sources in our survey. We 
have adopted the 90% completeness limit in this paper, which corresponds to the following 
limiting magnitudes: mi7'iiovi/=21.5 mag, mi?i60Ty=20.5 mag, and m^207M=18.0 mag. 



3. Results 

We present the basic results from the photometry. Objects are placed in the HR diagram 
in Section 3.1 based on the effective temperature estimate derived through the ffux ratio 
between F165M and F160W and the dereddened J band magnitudes. In Section 3.2, we 
discuss the J-H versus J color-magnitude diagram and the H-K versus J-H color-color 
diagram. Throughout this section, we assume all but two very blue bright stars are cluster 
members. The two blue stars have J — H < 0.5 mag and are assumed to be foreground 
objects. We have no direct estimate of the number of field stars expected towards Mon R2 
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which would contaminate our sample at the sensitivity of these observations. As an indirect 
test we have downloaded the 2MASS point sources located in a 2 degree region centered on 
(1=214.3 b=-15.4), where there coordinates for Mon R2 are (1=213. 7,b=-12.6). The 2MASS 
H band luminosity function is a power-law down to the completeness limit of H~15 mag. 
Assuming the field star population would follow a power-law to faint magnitudes, we would 
expect ^ 45 field stars within our field of view down to H=20.5 mag. However, the dust 
in the molecular cloud acts as a screen limiting the number of background stars we will 
detect. Carpenter et al. (1997) found the average extinction from the molecular cloud to be 
Av = 50 mag in the central 45" region and Ay = 33 mag on average across the central 3' 
square region. If the field star luminosity function extrapolated from 2MASS is reddened by 
Ay = 30 mag, the surface density of field stars down to a limiting magnitude of H=20.5 mag 
is expected to be 1-2 objects per square arcminute. We have also estimated the field star 
contamination using the synthetic Galactic model by Robin et al. (2003). The molecular 
cloud has been placed at 830pc and assumed to extinct background sources by Ay = 30 
mag. The model predict that we should observe 10464 sources per square degree down to 
H=20.5 mag, where 6022 sources are predicted to be foreground sources. Thus, the model 
predicts 1.7 foreground objects and 1.2 background objects; the former consistent with our 
identification of 2 blue foreground objects, the latter in agreement with the extrapolation 
from the 2MASS data. 



3.1. Effective temperatures and luminosities for objects in Mon R2 

A water vapor absorption feature is present in the H band in the spectra of stars and 
brown dwarfs with temperatures at or below 4000 K (see e.g. AUard et al. (2000)). This 
feature is covered by the F160W filter and since the F165M filter measures the continuum, 
the color mpiebM — "m-Fimw for an object cooler than ~ 4000 K can be used to estimate the 
effective temperature. We have compared the mpiQ^M — fnFiaow colors, with the synthetic 
colors predicted by the atmosphere models of Hauschildt et al. (1999) and Allard et al. 
(2000) for objects above and below 2500 K, respectively The mFi&bM — ^Fimw versus 
^FiiovK — 'mpimw color-color diagram is shown in Fig. 2. 
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Table 3: The fraction of recovered artificial stars as a function of magnitude for each filter. 
The detection limit was set to 5a. 



Mag 


FllOW 


F160W 


F165M 


F207M 


17.0 








97% 


17.5 








97% 


18.0 








96% 


18.5 








77% 


19.0 








9% 


19.5 




98% 


96% 


1% 


20.0 


99% 


98% 


98% 




20.5 


99% 


94% 


97% 




21.0 


97% 


44% 


51% 




21.5 


87% 


5% 


2% 




22.0 


66% 


2% 


1% 




22.5 


10% 








23.0 


2% 
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n^FiesM^nHpieow 

Fig. 2. — The 'mFi65M — "itt^figow versus uifiiow —'itT'Fwow color-color diagram. Overplotted 
are the DUSTY models by AUard et al. (2000) for temperatures below 2500 K (dashed hne) 
and the NEXTGEN models by Hauschildt ct al. (1999) for temperatures above 2500 K (solid 
line). The arrow indicates the magnitude and slope of Av=5 mag. The dotted lines indicate 
the region of the color-color diagram where a reliable temperature can be obtained if the 
photometric error in mFi65M —mpwow is less than ~0.06 mag. Objects for which the effective 
temperature has been derived are shown as filled circles and the error bars are indicated for 
these objects. The effective temperatures at positions on the two models are shown in units 
of 100 K and marked on the model with open diamonds. 
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The reddening vector is parallel to the models for temperatures above ~3800 K (mpiesM— 
TTT-Fwow > -0-1 mag) and no reliable temperature can be derived. There is a plateau from 
~3800 K and to ~2700 K (spectral types MO to M7.5) where an unique temperature can be 
obtained by projecting an observed point along the reddening vector until it intersects the 
temperature locus. The reddening vector crosses the models at two places for lower tem- 
peratures and a unique temperature cannot be obtained. We have estimated the effective 
temperature for the sources located between ~3300 K and ~2700 K which have magnitude 
errors smaller than 0.06 mag. An uncertainty of 0.06 mag translates into an uncertainty of 
the temperature determined from the model of ~10% for this temperature range. We have 
adopted the models with a surface gravity of logg — 4, consistent with typical values of 
low mass PMS objects. Wc do not necessarily trust the magnitude of the extinction needed 
for the observations to intersect the temperature locus, due to uncertainties in the pseudo 
continua predicted by the models. Previous studies have shown the strength of the water 
bands in the J and K bands only depend weakly on the surface gravity (Wilking et al. 2004; 
Gorlova ct al. 2003). It could thus be expected the H band water vapor absorption feature is 
also weakly dependent of gravity. Indeed, adopting a model with logg = 4.5 only decreases 
the derived temperature by < 50 K. Decreasing the surface gravity to logg — 3.5 has a more 
complicated effect. Whereas the low temperature regime (< 3000 K) is unaffected, we would 
underestimate the temperature by ~300 K by incorrectly adopting the logg = 4.0 model 
for objects with temperatures between 3000 K and 3500 K. Leggett et al. (1996) compared 
the NextGen models to observed spectra of late type stars. They found the predicted water 
bands to be deeper than observed for a given spectral type. The effect is that the models 
will predict a slightly too high temperature for a given waterband strength. Najita et al. 
(2000) and Gorlova et al. (2003) find reasonable agreement between temperatures derived 
from comparison of water vapor indices from the models and spectral types. However, Luh- 
man et al. (2003) found a systematic shift of spectral types relative to Najita et al. (2000) of 
1-2 subtypes which they attributed to difference in the surface gravity between field dwarfs 
used to calibrate the index of pre-main sequence objects. Based on the model prediction 
of the temperature accuracy and the indirect evidence from other water bands, we estimate 
our systematic uncertainties to be ±300 K for temperatures derived using the water vapor 
index. 

Objects with a derived effective temperature in the range 2700-3300 K and with a 
photometric error mpiQ^M — tT^Fimw < 0.06 mag are marked with filled circles in Fig. 2. We 
determined the bolometric luminosity adopting formula 2 from Gorlova et al. (2003), which 
assumes Mboi = 4.64, and BCv© = —0.19. The dominant source of error in the bolometric 
luminosity stems from uncertainties in the derived extinction. An error per filter of 0.06 
mag will translate into an error in the visual extinction Ay ~0.75 mag. Taking into account 
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an error of ~10% in the bolometric correction, and the distance uncertainty, we estimate 
the error in the logarithm of the bolometric luminosity to be 0.2. The objects are then 
placed in the HR diagram, which is shown in Fig. 3. We have not plotted source 83 in 
the HR diagram. Source 83 is a probable foreground source due to its very low extinction. 
Source 6 also appears under-luminous. One possible reason for the under-luminosity can be 
the extinction has been underestimated due to unresolved scattered light (cf. Wilking et al. 
2004). 

The objects from Carpenter ct al. (1997) with effective temperature and luminosity 
estimates arc shown as well. The 1, 10, and 1000 Myr Baraffc ct al. (1998) isochrones are 
shown. Since the Baraffe et al. (1998) models only extends to 1.4 Mq we have overplotted 
1 Myr and 10 Myr Palla & Stabler (1999) isochrones to cover masses up to 7 Mq. Table 4 
summarizes the parameters for the sources with the effective temperature derived from the 
water band method described above. 
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Fig. 3. — The HR diagram based on objects where the effective temperature and luminosity 

have been derived in this study using the mpiQ^M ~ ^Fimw color (open circles) and ob- 
jects from Carpenter et al. (1997) (filled circles). Overplotted are the 1, 10, and 1000 Myr 
isochrones made from BarafTe et al. (1998) (solid lines) and both Palla & Stabler (1999) 1 
and 10 Myr isochrones (dashed lines). 1 M0, 0.1 M©, and 0.04 M0 have been marked on 
the Baraffe et al. isochrones with plus signs. The Palla & Stabler (1999) isochrone covers 
the mass range 0.1-7 Mq. A typical error bar is shown. 
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Although a large fraction of the objects below 1 M© are located close to the 1 Myr 
isochrone, there appears to be a large scatter. Some the of scatter is due to errors in 
estimates of the effective temperatures and luminosities and some of the scatter might be 
real. However, most of the objects that appear to be older than 1 Myr are more massive 
than 1 Mq, which is the maximum mass we attempt to constrain the IMF. We find that 
for the objects in the HR diagram with masses between 0.04 Mq and 0.4 M0 that 7 objects 
are 1 Myr or younger and 6 are older than 1 Myr. A median age of 1 Myr appears to be 
appropriate for the lower mass content in Mon R2. 
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Table 4: Physical parameters for the sources where we have estimated the effective temper- 
ature. 



ID 


Av" 


Tcff 


Mj 


excess^ 


logL 


3 


10.4 


3250 


5.29 


N 


-1.00 


6 


1.4 


2950 


9.35 


Y 


-2.66 


13 


9.4 


3250 


5.59 


N 


-1.12 


43 


21.1 


2750 


2.93 


Y 


-0.14 


82 


4.6 


2950 


4.00 


Y 


-0.51 


83 


0.0 


3150 


10.54 


Y 


-3.11 


105 


17.7 


3150 


6.18 


Y 


-1.37 


144 


20.4 


2950 


6.81 


Y 


-1.64 


147 


14.2 


3150 


3.90 


N 


-0.45 


151 


16.7 


3250 


5.37 


Y 


-1.03 


165 


9.7 


3150 


3.33 


Y 


-0.28 


171 


6.5 


3250 


5.26 


N 


-1.00 



"Derived from the color-magnitude and color-color diagrams in Figs. 4 and 5 as described in Section 3.2. 
''Estimated from the color-color diagram in Fig. 5. 
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3.2. Color-magnitude and color-color diagrams 

We present in Fig. 4 the J-H versus J color-magnitude diagram for the the inner 1' 
square of Mon R2. 
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Fig. 4. — The J-H versus J color-magnitude diagram for Mon R2 based on data from Table 2 
converted into the CIT system. Overplotted as the solid lines are a 1 Myr isochrone from 
the Baraffe et al. (1998) models as discussed in the text, and the same isochrone reddened 
by Ay = 19 mag. The symbol sizes indicate photometric error, where the error in the 
transformation into the CIT system has been included. The horizontal line indicates the 
90% completeness hmit of the observations. The arrow illustrate the effect Ay = 5 mag 
extinction. Dashed lines are drawn between the unreddened and reddened isochrones at 
object masses of 1, 0.1, and 0.08 Mq. The location of an unreddened 0.02 Mq object is 
marked. 
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Also shown are a one Myr isochrone created from the Baraffe et al. (1998) models, both 
unreddened and reddened by Ay = 19 mag and shifted with a distance modulus of 9.6 mag 
(830 pc). The Baraffe et al. models with Lmix = Hp are chosen since they cover the whole 
mass range where we attempt to constrain the IMF (0.02 Mq-I.O Mq) and they are tied to 
the atmosphere models used in Section 3.1. 

In order to convert the effective temperatures and bolometric magnitudes provided by 
the models to observables we have used the main sequence colors from Bessell & Brett (1988) 
and the temperature scales from Schmidt-Kaler (1982) and Bessell (1991) for objects earlier 
and later than spectral type than K7, respectively. Data from Dahn et al. (2002) have 
been used for spectral types later than M6. The bolometric corrections and colors were 
interpolated using spline interpolation. Due to large scatter observed in the intrinsic colors 
and magnitudes for the late-type objects, a linear fit has been performed for objects later 
than M6. 

We expect the observed color of pre-main sequence objects to be due to several effects. 
One is the intrinsic color of the object that depends on mass and age, which in turn fix Teff 
and log(g'). Also, the general interstellar extinction towards the object and the possible pres- 
ence of a circumstellar disk associated with the object. Prom the color-magnitude diagram 
it is not possible to distinguish between the two latter effects. However, they can, at least 
partly, be separated in the H-K versus J-H color-color diagram. Objects without strong 
emission from a disk have colors resembling reddened main sequence stars and therefore 
populate a confined region in the color-color diagram. Fig. 5 shows the color-color diagram 
for Mon R2 in the CIT system. Overplotted are the colors for main sequence stars from 
Bessell & Brett (1988), reddening vectors from stars with spectral type M6 and K7, and the 
dereddened classical T-Tauri locus derived by Meyer et al. (1997). Objects with an excess 
have a larger reddening in H-K relative to what would be expected from their J-H color and 
are located to the right of the main sequence stars in the color-color diagram. 
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Fig. 5. — The H-K versus J-H color-color diagram for Mon R2 based on data in Table 2 
converted into the CIT system. The open circles indicate stars not identified with any 
near-infrared excess whereas the filled circles indicate stars with near-infrared excess and an 
estimated spectral type earlier than M6 (see text). Shown as the dotted line are the colors 
from Bessell & Brett (1988) for main sequence stars with spectral types B8 to M6. The 
length of the arrow illustrates the effect of Ay = 5 mag extinction. The dashed line is the 
unreddened T-Tauri locus from Meyer et al. (1997). Extended from the main sequence colors 
and the classical T-Tauri slope are solid lines parallel to the reddening vector. 
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A large fraction of the stars fall in the region expected for stars with no infrared excess 
and spectral type earlier than M6. The objects having larger than expected H-K relative 
to the reddening vector extending from the M6 main sequence star can either be stars with 
spectral type earlier than M6 with an infrared excess or they can be objects of later spectral 
type with or without infrared excess. The intrinsic colors for very late type objects (later 
than M6) almost overlap in the color-color diagram with the dereddened classical T-Tauri 
locus. Thus, we cannot tell for late spectral type objects whether they have a disk or not 
from the color-color diagram alone. 

The error adopted in J — by not being able to distinguish between an excess object 
and a late type object will be < 0.1 mag, corresponding to Ay < 1 mag. We have therefore 
dereddened all objects in this part of the diagram to the T-Tauri locus in order to estabhsh 
the intrinsic magnitudes of the objects. For stars with colors consistent with main sequence 
stars they can in general be deredden to two different intersections with the main sequence 
colors. Therefore, we have dereddened these stars using the J-H color-magnitude diagram 
adopting the 1 Myr isochrone (Meyer 1996; Carpenter et al. 1997; Wilking et al. 2004). 

4. Analysis 

In this section, we use the HST photometry to constrain the stellar and sub-stellar 
IMF in the central part of the Mon R2 cluster as well as re-visit the fraction of stellar mass 
objects with circumstellar disks. We begin by first defining an extinction limited sample 
that is complete down to 0.1 Mq. We use this sample to determine the fraction of stars 
with circumstellar disks and to search for variations in the extinction as a function of object 
mass. We then define three extinction limited samples complete for both stars and brown 
dwarfs in order to constrain the sub-stellar IMF of Mon R2 and compare these ratios with 
the field star IMF and other star forming regions. 

4.1. The circumstellcir disk fraction for stcirs in Mon R2 

What is the fraction of stellar objects in our extinction limited sample showing evidence 
for a circumstellar disk as discussed above? We have constructed an extinction limited 
sample in order to obtain a representative sample across the whole mass range considered. 
Without an extinction limited sample, we would underestimate the number of fainter, lower 
mass sources relative to intrinsically more luminous objects of higher mass. The maximum 
extinction is determined such that the sample considered will be complete down to spectral 
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type M6, corresponding to a 0.1 star for the 1 Myr isochrone. This corresponds to an 
extinction hmit of Ay = 19 mag. We have excluded the two objects on the bhic side of the 
isochrone in Fig. 4 since these objects are probable foreground stars (star ID 61 and ID 83). 
An exclusion of these objects is equivalent to excluding objects with a derived extinction 
less than Ay = 1 mag. The extinction limited sample (Ay = 19 mag) contains 43 objects 
in total. The selected sources are all located between the two upper dashed lines in Fig. 4. 
In Section 4.4 we construct three further extinction limited samples in order to address the 
IMF into the brown dwarf regime. 

For objects with a mass between 0.1 M0 and 1.0 Mq, wc find a disk fraction of 27±9%. 
For a wider area but more shallow survey of Mon R2, Carpenter et al. (1997) found a disk 
fraction of 48±8% from JHK photometry. A disk fraction of ~ 30% is relatively low compared 
to what is found in other star clusters of comparable age. For example, in the Orion Nebula 
Cluster (ONC), Lada et al. (2000) found a disk fraction of 50 ±20% using JHK photometry, 
(cf. Hillenbrand et al. 1998), which within the error bars is consistent with Mon R2. Haisch 
et al. (2001) found a disk fraction of 21 ± 5% in IC 348 using JHK photometry. 

4.2. Constraints on the IMF in the center of Mon R2 

With a dataset significantly deeper than the previous study by Carpenter et al. (1997), 
we can constrain the IMF into the brown dwarf regime. The use of an extinction limited 
sample for IMF studies assumes the extinction distribution is independent of the objects 
mass. If this is not the case, a bias would be introduced by preferentially excluding either 
high or low mass objects. To test if the extinction distribution is independent of object 
mass, we have divided our sample into objects between 0.2-1 Mq and between 0.03-0.2 
Mq for an Ay < 10 mag extinction limited sample (see below). We have then performed a 
two sided Kolmogorov-Smirnov (KS) test whether the extinction distributions are consistent 
with having been drawn from the same parent population. The KS test returned a value of 
0.49, from which we conclude there is no obvious difference in the extinction distribution as 
a function of mass. 

We have calculated the ratio of stars between 0.08 and 1 Mq divided by the number of 
objects between 0.04 and 0.08 Mq, i?=N(0.08-1.0 Mo)/N(0.04-0.08 Mq), for an extinction 
fimited sample Ay < 13 mag based on 34 objects. The ratio is found tobei?=10.3±5.8 
for a cluster age of 1 Myr, where the errors are derived assuming Poisson errors. We then 
compared the ratio with the similar predicted ratio from the Chabrier (2003) system field 
star IMF (dN/dlogMoc exp {-{login - log{).22f/{2 ■ 0.57^))). The distribution of predicted 
ratios is given by the binomial distribution and is shown in Fig. 6. The height of each line 
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Fig. 6. — A comparison of our derived R value with the R values obtained from a Chabrier 
(2003) field IMF for a sample size of 34 objects. The predicted ratio distribution from the 
field IMF was determined used the binomial distribution. 
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Although the ratio derived for Mon R2 is slightly higher than the peak of the distribution 
of ratios predicted by the Chabrier (2003) system IMF, the probability of obtaining the 
derived value or higher from the Chabrier system IMF is 19%. 

We have constructed two additional extinction limited samples, one down to Ay — 10 
mag and one to Ay = 7 mag. The two samples are complete down to 30 Mjup, and 20 
Mjup. respectively. The upper mass for both samples is still 1 Mq. Although both samples 
extend to fainter absolute magnitudes than the Ay = 13 mag sample, the lower maximum 
extinction reduce the sample sizes significantly. Only 19 objects are included in the Ay = 10 
mag sample and 13 in the Ay = 7 mag sample. Wc have calculated the ratio of high- 
mass to low-mass objects for the two latter samples in a similar manner as above. We find 
the ratios i?*=N(O.O8-l.OM0)/N(O.O3-O.O8M0)=8.5±6.4, and i?**=N(O.O8-l.OM0)/N(O.O2- 
0.08Mq)=2.2 ± 1.3. The similar expected ratios from the Chabrier (2003) system IMF are 
R* — 4.2, and R** — 3.5, respectively. It appears that the ratio of stars to brown dwarfs 
observed relative to the ratio expected from Chabrier (2003) decreases (/?=10.3±5.8 vs 
5.3;R* = 8.5 ± 6.4 vs 4.3; R** = 2.2 ± 1.3 vs 3.5) as a function of depth into the cloud. If 
confirmed, this might suggest that brown dwarfs and stars are not uniformly distributed in 
MonR2. However, the KS test did not indicate the distribution of extinction is different for 
the high and low mass samples, down to 0.03 Mq, so we do not believe these differences 
to be significant. All the observed ratios are consistent with the Chabrier (2003) IMF. The 
cumulative distributions of the absolute J band magnitudes for the three extinction limited 
samples are shown in Fig. 7. The location of a 0.08 Mq 1 Myr object from the Baraffe et al. 
(1998) models is indicated. 
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Fig. 7. — The cumulative distributions of the absolute J band magnitudes for the objects 
in the three extinction limited samples, Ay = 13 mag (sohd hne, 34 objects), Ay = 10 mag 
(dotted hne, 19 objects), and Ay = 7 mag (dashed hne, 13 objects). The absolute magnitude 
of a 1 Myr 0.08 Mq object (Baraffe et al. 1998) is marked with the vertical dot-dot-dashed 
line. 
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5. Discussion 

The derived ratios of high-mass to low-mass objects can be compared with similar 
ratios derived from other young star clusters with similar spatial resolution. We will here 
focus on some of the best studied clusters, namely the ONC, IC348, and Taurus and compare 
the values of R, R*, and R** for those clusters with the values derived here. We will briefly 
discuss the selection criteria for the IMF samples for the different clusters and we compile 
the ratios in Table 5. 

As the closest site of massive star formation, the ONC has been observed in some 
detail. The IMF has been estimated through deep near-infrared imaging (Hillenbrand & 
Carpenter 2000; Muench et al. 2002). Using H and K band observations, Hillenbrand & 
Carpenter (2000) derived the IMF down to 0.02 Mq for Ay < 10 mag. From their histogram 
we have derived the ratios for the ONC where we find R = 5.5 ± 0.8, R* = 4.3 ± 0.6, and 
R** = 3.0 ±0.6, respectively. In another imaging study, the IMF of the ONC was constrained 
through model fits to the K band luminosity function by Muench et al. (2002) down to ~10 
Mjup. The ratio of stars to brown dwarfs was found by integrating the functional fit to the 
IMF presented by Muench et al. (2002). From the functional fit, we find R = 6.4, R* = 4.7, 
and R** — 3.7, respectively consistent both with result by Hillenbrand & Carpenter (2000) 
and our ratios for Mon R2. 

Spectroscopic observations of objects in the ONC by Slesnick et al. (2004) indicate the 
ratio of high mass to low mass objects might be higher than found from imaging alone. 
A ratio of R'=N(0.08-0.4 Mq)/N(0.03-0.08 M©) =4.5±0.7 is found from the spectroscopic 
survey. It should be kept in mind that correction factors of ~2 have been applied to the 
number counts as discussed in Slesnick et al. (2004). The comparable ratio for Mon R2 is 
4.5±3.5 which is in excellent agreement with the Slesnick et al. (2004) result. For comparison, 
the similar ratio from the Muench et al. (2002) IMF is 4.1. 

As one of the nearest low-mass star forming regions to the Sun, Taurus has been inves- 
tigated in detail. Luhman (2004) has presented the IMF down to 0.015 Mq for an extinction 
limited sample Ay < 4 mag. We have calculated the ratio of high to low-mass stars from 
the data presented by Luhman (2004). We find the derived IMFs for Taurus and Mon R2 
to be consistent within the uncertainties for all three ratios. 

The ratios for Taurus might be suspect for two reasons. Luhman (2004) emphasized 
the sample is incomplete for the mass interval 0.3-0.6 Mq due to saturation, resulting in the 
derived ratios being lower limits. Also, the surveyed region by Luhman (2004) only covered 
parts of Taurus (4 square degrees), preferentially centered on the dense filaments. A larger 
area survey (28 square degrees) has found relatively more brown dwarfs (Guieu et al. 2006), 
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although the latter survey only extended down to 0.03 Mq, for an extinction limited sample 
Ay < 4 mag. They concluded the ratio of brown dwarfs to star in Taurus and the ONC 
arc similar within the statistical uncertainties. Although the ratio derived by Guicu ct al. 
(2006) extended to higher masses relative to this paper, we can still estimate the effect of 
their updated number of brown dwarfs in Taurus. The newly detected brown dwarfs in the 
larger surveyed region by Guieu et al. (2006) reduced the derived ratios by ~ 25% relative 
to Luhman (2004). Correcting the Taurus results quoted in Table 5 downwards by 25% only 
improves the agreement with the ratios for Mon R2. 

The low-mass stellar content of 1C348 has been investigated recently by Luhman ct al. 
(2003). The IMF was derived down to 30 Mjup, using an extinction limited sample Ay < 4 
mag. The objects presented by Luhman et al. (2003) have been spectroscopically confirmed 
as members. The spatial resolution in the imaging survey was O'.'G. Scaled to the distance of 
Mon R2, this is a similar physical resolution as our F207M band image. Based on the data 
published in Luhman (2004), we have derived the ratios R = 16.8 ±5.8 and R* = 11.6 ±3.3. 
The ratios for IC 348 appears to be relatively high but still consistent within 2o" to the 
results for Mon R2. The IMF has been probed in IC 348 down to 0.015 M0 utilizing 
HST/NICMOS imaging (Najita et al. 2000). Using the the strength of the 1.9/im water 
band they determined the spectral type of objects in the mass interval 0.015-0.7 Mq. We 
have calculated the modified ratio of high-mass to low-mass objects using objects between 
0.25-0.7 M0 as the high-mass objects and 0.02-0.25 Mq as the low-mass objects. The ratio 
derived from the histogram in Najita et al. (2000) is 0.4, whereas the similar ratio for Mon 
R2 is 0.3±0.2 for an extinction limited sample. Ay < 7 mag. 

One issue that should be explicitly taken into account in these comparisons is the role of 
binaries. It is well known that a large fraction of stars reside in multiple systems and due to 
our finite angular resolution, we can only resolve a certain fraction of the binary systems in 
each cluster. The resolution in the F207M band with HST/NICMOS at a distance of 830 pc, 
resolves all systems more widely separated than ~160 AU. Both the ONC imaging surveys 
have a similar physical resolution with the cluster at 430 pc. Luhman (2004) considered all 
objects with a separation smaller than 2" in Taurus as single objects. For a distance of 140 pc, 
this corresponds to a separation of 280 AU. On the other hand, the HST/NICMOS imaging of 
IC 348 by Najita et al. (2000) has twice the physical resolution as this study, resolving binaries 
with separation larger than ~ 75 AU. Since the physical angular resolution is similar within 
a factor of ~2 for most of the clusters compared, we do not expect unresolved binaries to 
effect the ratios differently for the clusters discussed here. However, the field IMF by Kroupa 
(2002) has been constructed from individual objects. The majority of binaries in young open 
clusters have separations smaller than ~ 100 AU (Patience et al. 2002, cf. Duquennoy & 
Mayor (1991)) and would be unresolvable in our observations. The quantitative change of 



-32- 



the IMF due to unresolved binaries depends on the details of the binary frequency versus 
separation and the distribution of mass ratios, neither of which are well known, especially 
for low masses. However, there is some indication the binary frequency is lower for cool 
objects (Burgasser et al. 2006). Chabrier (2003) have presented the system IMF assuming 
no unresolved binaries derived from the local 8 pc solar neighborhood field star sample, 
which is used in Table 5. 

The derived ratios for Mon R2 appears consistent with the similar ratios for the ONC, 
IC 348, and Taurus. They are also consistent with the ratios derived for both the Kroupa 
(2002) and Chabrier (2003) IMF. Thus, despite more than an order of magnitude difference 
in total cluster mass, it appears the IMF down to 0.02 M© is similar. So far, there is little 
evidence for variations in the sub-stellar IMF at least down to 20 Mj^p. It remains to be 
seen whether or not the ensemble of the observations can distinguish between a flat or a 
faUing IMF between 0.02-0.08 M© (cf. Allen et al. 2005). 

Several theoretical considerations would predict a deflcit of brown dwarfs in Taurus rel- 
ative to more massive regions. Goodwin et al. (2004) suggest the IMF in a region like Taurus 
should have relatively fewer brown dwarfs due to the narrow distribution of core masses in 
Taurus relative to regions hke the ONC. An alternative explanation has been presented by 
Bate & Bonnell (2005). They suggested the lower density in the Taurus molecular cores 
results in a higher Jeans mass which in turns results in a higher average mass. Goodwin 
et al. (2004) indicate through numerical simulations that the peak of the IMF will shift to 
lower masses as the degree of turbulence increases. 

The increase of the relative number of brown dwarfs when the minimum mass of the 
sample is decreased is curious. Since the extinction limit decreases as the limiting mass is 
extended to lower masses, one possibility is that the more massive stars are observed prefer- 
entially deeper within the molecular cloud. However, we found no evidence for variation in 
the extinction as a function of object mass down to 0.03 M©. One possibility, although spec- 
ulative, is that the brown dwarfs have been ejected from small N body systems as proposed 
by e.g. Reipurth & Clarke (2001). We would then detect the brown dwarfs ejected towards 
us (preferentially with lower extinction) even though the parent system containing the star 
would be located deeper in the cloud. Since we are only probing the surface of the cluster 
in the most shallow extinction limited samples, such a scenario would explain the observed 
trend. Deeper analysis of the dynamical evolution (velocities and spatial distribution) of 
cluster members as a function of mass is needed. 
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6. Conclusions 

We have presented the results from HST/NICM0S2 FllOW, F160W, F165M, and 
F207M band imaging of the inner 1' x 1' of the embedded cluster associated with Mon 
R2. Our results are as follows: 

• The effective temperature has been estimated for a small set of stars in the temperature 
range 2700 — 3300 K based on a water vapor index and these objects have been placed 
in the HR diagram. A 1 Myr Baraffe et al. (1998) isochrone is consistent with the lower 
mass objects placed in the HR diagram in agreement with Carpenter et al. (1997). 

• We find that for stars with spectral type M6 or earher (O.l-IM©) a disk fraction of 
27 ± 9% based on an extinction hmited sample of 43 stars Ay < 19 mag. 

• We created three extinction limited samples complete for Ay = 13, 10, and 7 mag, 
and containing 34, 19, and 13 objects, respectively. We have calculated the three 
ratios of low mass stars to brown dwarfs R = iV(0.08 - l.OM0)/(O.O4 - O.OSM©) = 
10.3 ± 5.8, R* = N{0.08 - l.OM0)/iV(O.O3 - O.O8M0) = 8.5 ± 6.4 and R** = A^(0.08 - 
l.OM0)/A^(O.O2 - O.O8M0) = 2.2 ± 1.3. 

• The derived ratios are consistent with the similar ratios for Taurus, IC 348, the ONC, 
and the system field IMF of Chabrier (2003). There is thus no compelling evidence for 
variations in the relative brown dwarf content between Mon R2 and other nearby star 
forming regions. 
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Table 5: The derived ratios of high-mass to low-mass objects for Mon R2, the field star IMF 

and ot her star forming regions. 

Mio^ MonR2 Field^ Field^ IC348" Taurus'^ ONC" ONC^ 
~R OM 10.3±5.8 2^9 5^3 16.8±5.8 9.6±3.2 6A 5.5±0.8 
R* 0.03 8.5±6.4 2.2 4.2 11.6±3.4 6.9±2.0 4.7 4.3±0.6 
R** 0.02 2.2±1.3 1.7 3.5 5.4±1.4 3.7 3.0±0.6 



"Kroupa (2002), single star IMF 
''Chabrier (2003), system IMF 
'^Luhman et al. (2003) 
•^Luhman (2004) 
^Muench et al. (2002) 
■''Hillenbrand & Carpenter (2000) 



